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ABSTRACT 

Context. Jets from the youngest protostars are often detected only at mm wavelengths, through line emission of CO and SiO. However, 
it is not yet clear if such jets are mostly molecular or atomic, nor if they trace ejected gas or an entrained layer around an embedded 
atomic jet. 

Aims. We investigate the warm gas content in the HH21 1 protostellar outflow to assess the jet mass-flux in the form of H 2 and probe 
for the existence of an embedded atomic jet. 

Methods. We employ archival Spitzer slit-scan observations of the HH21 1 outflow over 5.2-37 //m obtained with the low resolution 
IRS modules. Detected molecular and atomic lines are interpreted by means of emission line diagnostics and an existing grid of 
molecular shock models. The physical properties of the warm gas are compared against other molecular jet tracers and to the results 
of a similar study towards the L1448-C outflow. 

Results. We have detected and mapped the v=0-0 S(0) - S(7) H 2 lines as well as fine-structure lines of S, Fe + , and Si + . The H 2 is 
detected down to 5" from the source and is characterized by a "cool" T~ 300K and a "warm" T~ 1000 ± 300 K component, with an 
extinctionAv ~ 8 mag. The amount of cool H 2 towards the jet agrees with that estimated from CO assuming fully molecular gas. The 
warm component is well fitted by C-type shocks with a low beam filling factor ~ 0.01-0.04 and a mass-flux similar to the cool H 2 . 
The fine- structure line emission arises from dense gas with ionization fraction ~ 0.5 - 5 x 10~ 3 , suggestive of dissociative shocks. 
Line ratios to sulfur indicate that iron and silicon are depleted compared to solar abundances by a factor ~ 10-50. 
Conclusions. Spitzer spectral mapping observations reveal for the first time a cool H 2 component towards the CO jet of HH211 
consistent with the CO material being fully molecular and warm at ^ 300 K. The maps also reveal for the first time the existence of 
an embedded atomic jet in the HH21 1 outflow that can be traced down to the central source position. Its significant iron and silicon 
depletion excludes an origin from within the dust sublimation zone around the protostar. The momentum-flux seems insufficient to 
entrain the CO jet, although current uncertainties on jet speed and shock conditions are too large for a definite conclusion. 

Key words. STARS: FORMATION, ISM: JETS AND OUTFLOWS, ISM: INDIVIDUAL OBJECTS: HH211-mm, INFRARED: 
ISM: LINES AND BANDS 



1. Introduction 

The process of mass accretion leading to the formation of low 
mass protostars is always associated with ejection of material 
in the form of well collimated high-velocity jets and/or wider 
slower bipolar outflows (Rei purth & Bachiller|1997)|Arce et"aT] 
|2007| l. Accretion and ejection phenomena are believed to be in- 
timately connected through the existence of a magnetized accre- 
tion disk around the forming protostar, but there is no consensus 
so far about the exact ejection mechanism (Pud ritz et al.|2 007 



Shang et al. 2007 Ferreira et al. 2006). The study of the physica 
properties of jets (e.g. temperature, density, abundances, mass 
flux) is therefore crucial to obtain indirect constraints on the 
mass accretion and ejection processes, especially in the case of 
protostars in the earliest evolutionary stages - the so called Class 



protostars - where the launching zone is heavily embedded in 
a dense cocoon of dust and gas. 

Due to this high extinction, jets from the youngest Class 
sources cannot be seen in the optical and are mostly traced in 
near-IR H 2 and in mm CO and SiO lines ([Guilloteau et al.|1992[ 



McCaughrean et al.|1994||Davis & Smith|1996[ ). However, it is 
still unclear whether the molecules trace a sheath of entrained 
ambient material around an unseen underlying atomic jet, or if 
they trace the primary jet material itself. In the latter case, the 
chemical composition of the jet (molecular vs. atomic fraction, 
dust depletion) holds crucial clues to the ejection zone, as the 
strong UV flux generated by the accretion shock should destroy 
dust grains and H2 molecules in the innermost regions of the 
disk. A wind from such inner zones would be mostly atomic even 
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if CO and SiO are abundant — unless the mass flux is unusually 
large ( |Ruden etal.|1990) . 

Mid-infrared spectroscopic observations are essential to 
progress on these questions. A Spitzer study of the innermost 
jet regions of the Class L 1448 outflow ( |Dionatos et al.|2009| > 
revealed an underlying deeply embedded atomic/ionic compo- 
nent, not seen in the near-IR range, as well as mid-IR H2 emis- 
sion tracing the warm molecular content. Here we present a more 
extensive Spitzer mapping study of the HH21 1 outflow, a partic- 
ularly interesting target that appears as a "text-book" example of 
jet-driven flow. 

HH211 is located in the IC 348 complex in Perseus at a 
distance estimated to be between 250 pc (Eno ch et al.| 2006), 
adopted in the current paper, and ~320 pc (Herbig 1998; Lada 
et al.|2006) . The driving source HH211-mm is detected in mm 



continuum emission (e.g. Gueth & Guilloteau 1999; Lee et al. 



2007 ), and has been classified as a low-mass and low-luminosity 
Class young stellar object ( Froebrich 200 5} . 



The outflow was discovered by McCaughrean et al. ( 1994 ) 
in near-IR H2 emission tracing 2000 K hot, shocked gas in two 
bright symmetric bowshocks separated by 0.1 3pc with a faint 
chain of knots in between. A few compact knots of atomic emis- 
sion at optical and near-IR wavelengths (i.e. H a , [SII]6730A 
and [Fell] 1 .64//m) are also seen near the brightest H2 peaks 
O'Conne ll et al. 12005 ||Caratti o Garatti et al.|2006HWalawender| 



let aLJ2005||2006|l. 



The absence of any further shock structures on larger scales 
makes it the youngest outflow known so far, with a kinematical 
age of only 1000 yr x(V /100) km s _1 ( jEisloffel et aI1|2003a ' 
IGueth & Guilloteau|1999| l. 



Lo w velocity CO J= 2-l ( |Gueth & Guillotea u[T999| and CO 



J=3-2 (Lee et al. 2007 1 observations delineate the shape of a 



bipolar swept-up cavity with its ends tracing the H2 bow-shocks 
and its flanks connecting back to the driving source, while the 
high- velocity components of the same transitions reveal an in- 
ner well-collimated CO jet extending out to +25" from the 
source. The CO jet is also traced in SiO through interferomet- 
ric observations of the J= 1-0 ([Chandler & Richer|[200T] l, 5-4 
( |Hirano et aLl|2006j l and 8-7 ( [Lee et al.||2007|l transitions, and 
single-dish multi-transition studies (Gibb et al. 2004; Nisini et al. 



2002b I. Excitation of such high J levels implies that the SiO 



jet is warmer and denser than the outflow cavity traced by low- 
velocity CO emission. 

All observations reveal a highly bipolar structure, with blue 
and red-shifted lobes pointing southeast and northwest from the 
driving source. The two lobes are well separated and such mor- 
pholog y combined with the low o bserved radial velocities (~20 
km s^ jGuefh & Guilloteau|l9 99) implies that the inclination of 
the outflow is less than 10° from the plane of the sky (Gueth & 



|Guilloteau|1999l|Chandler & Richer|2001l|Leeet al. 2007 ). This 
simple geometry is favorable to modelling. Indeed, the CO cav- 
ity shape is well fitted by a dynamical model of jet-driven bow- 
shock propagating into a medium of decreasing density (Gueth 



& Guilloteau| [1999). This interpretation is supported by fitting 



the near-IR line fluxes and morphologies at the apex of the red- 
shifted lobe with a series of bowshocks ( lO'Connell et al.|2005| l. 
A Spitzer map obtained at the tip the blue lobe indicates a higher 
excitation 40 km/s dissociative shock on that side, with a sub- 
stantial UV flux ( [Tappe et al.|2008] >. 

In this paper, we present archival Spitzer-IRS slit-scan ob- 
servations covering most of the HH2 1 1 flow with low spectral 
resolution (Section 3), which we analyse to address the physical 
conditions and the molecular and dust content. In Section 4, the 
detected rotational H2 transitions are employed as probes of the 




Fig. 1. Continuum subtracted spectrum, integrated over a region 
of 612 square arcseconds covering the HH211 flow (middle), 
shown in comparison with a spectrum of the off-outflow region 
of the same area (top). H2 emission from the 0-0 S(0)- S(7) 
lines is detected in the outflow, together with atomic emission 
from the ground state transitions of [Fell], [Sill] and [SI]. The 
H2 S(0) and S(l) lines are detected also in the off-outflow spec- 
trum, indicating the presence of a diffuse line emission compo- 
nent not related to the HH21 1 flow. Diffuse emission from PAH 
features, such as the 1 1 3pm feature, is also evidenced. (Bottom:) 
On source spectrum extracted from an area equal to the LL pixel 
scale (110.25 square arcseconds), after removing the zodiacal 
light and ambient contribution presenting the intrinsic contin- 
uum emission from the outflow source. 



physical conditions in the jet close to the driving source and the 
shocked gas further out, and compared with an existing grid of C 
and J-type shock models. Atomic lines present in our spectra are 
treated separately as probes for the existence of a deeply embed- 
ded atomic jet, and used to measure dust depletion. Mass-fluxes 
are inferred for the H2 and atomic components, and compared to 
those measured in CO and SiO, as well as to a similar study of 
the inner jet of L1448. Our conclusions are presented in Section 
5. 



2. Observations and data reduction 



Observations were obtained with the Spitzer satellite (Werner 



et al. 2004) and retrieved from the Spitzer Public Data Archive 
using the Leopard software. They were performed as part of the 
"Shock dissipation in Nearby Star Forming Regions" program 
conducted by J. Bally (P.I.). In these, the low resolution modules 
(R~60-130) of the Spitzer Infrared Spectrograph (IRS, Houck 
|et al.||20"04"| ) were used in slit-scan mode to cover an area of 
57" x 73" and 157" x 168" for the Short Low (SL) and Long 
Low (LL) modules respectively, centered on the Class source 
HH211 -mm (a/2ooo : 03h43m56.5s, £/2ooo : 

+32d00m51s). A per- 
pendicular scan step equal to the slit width of each module was 
used, with a total integration time of 171 min. The combination 
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Fig. 2. H2 v=0-0 S(2) - S(7) line intensity maps of the HH211 flow observed with the Spitzer SL module, reconstructed with a 
pixel scale of 3.5". Contours are from 10~ 12 W cirT 2 sr with an increment of 8 10~ 13 W cm~ 2 sr -1 for the para lines (even J), 
and from 2 1CT 12 W cirT 2 sr~' with an increment of 2.5 1CT 12 W cirT 2 sr~' for the ortho transitions (odd J). Crosses indicate the 
driving source position at O720oo:03h43m56.5s, <5/200o:+32d:OOm:51s. In the H2 S(5) map, points B1-B3 and Rl, R2 indicate peaks 
of emission, where Ay and ortho to para ratio are estimated and excitation analysis is carried out. 



of both IRS low resolution modules gives a complete wavelength 
coverage of 5.2 - 37.0 yum. 

Initial data processing was performed at the Spitzer Science 
Center using Version 15 of the processing pipeline. Spectral 
data-cubes were built using the CUBISM software package 
(Smith et al. 2007) and bad/rogue pixels were masked through 
visual inspection. 

Figure [TJ presents the extracted continuum-subtracted spec- 
trum, integrated over a region encompassing the HH211 flow. 
The full series of H2 pure rotational lines (S(0)-S(7)) were de- 
tected, along with atomic and ionic lines from the fundamental 
transitions of [SI], [Fell], and [Sill] at 25/im, 26/mi, and 35yum 
respectively. For comparison, Fig.[T]shows also a spectrum of a 
region encompassing the same area but in a direction perpendic- 
ular to the outflow. The H2 S(0) and S(l) lines are detected also 
in this off-outflow spectrum, indicating the presence of a diffuse 
line emission component not related to the HH21 1 flow. Diffuse 
emission from PAH features, such as the 1 13(xm feature, is also 
evidenced. In Figure [T] we also plot the on source spectrum for 
an area equal to the LL pixel size, after subtracting the contribu- 
tion from an off-source, free of line-emission position. This only 
removes zodiacal light and reveals the intrinsic mid-IR contin- 
uum from the central source, which is typical of low-luminosity 
Class sources. 

Subsequent analysis consisted in the construction of individ- 
ual line emission maps, using a home-built pipeline; in this, for 
each spatial pixel of the data-cube, the brightness of each spec- 
tral line of interest was calculated by Gaussian fitting after sub- 
tracting a local second order polynomial baseline. The resulting 
line intensity maps have a square pixel of side equal to the slit 
width of the IRS module, namely 3.5" and 10.5" for the SL and 
LL modules respectively, while the diffraction limit of the tele- 
scope is 2.4" and 6.0" at /l=10yum and 25yt/m respectively. The 



astrometric accuracy of the maps was found to be good within 
the limits imposed by the pixel size of each IRS module. 



3. Spectral line maps 

3.1. H 2 emission 

Figure 2]presents the emission line maps of the S(2) - S(7) ro- 
tational H2 lines observed with the SL - IRS modules with a 
3.5" sampling. Contours shape a characteristic bipolar outflow 
pattern, where the H2 emission is detected down to a projected 
angular distance ~ 5" from the driving source. Further down- 
wind, peaks of emission which can be attributed to shocked gas 
are observed. In the H2 S(5) map of Fig. 2]we label these peaks 
as B1-B3 in the southeast, blue-shifted lobe and R1-R2 in the 
northwest, red-shifted lobe, respectively. 

Figure [3] shows the H2 S(5) emission map overlaid on the 



1999> 



high- velocity CO J = 2— >1 map of Gueth & Guilloteau 
(upper panel) and on the SiO J = 8— >7 and CO J = 3— >2 maps 
of Lee et al. (|2007) (lower panels). The grayscale background 
shows the near-IR H2 v=l-0 S(l) images from McCaughrean 
|et al.| ( |1994] > or |Hirano et al.| ( |2006) . Overall, the spatial corre- 
spondence between mid-IR and near-IR H2 lines is quite good 
except towards the broadest parts of the outflow cavity, where 
midIR emission appears to trail behind the near-IR one. 

In the southeastern - blue lobe, the peak of mid-IR H2 emis- 
sion closest to the source (Bl) coincides with the 2.12/jm H2 
knot G in the nomenclat ure of |McCaug hrean et al. ( 1994]), the 
BII knot of the CO j et ([Gueth' & Guilloteau||1999| l and knots 
BK4-6 in the SiO jet ( |Lee et al.|2007| > (see lower panel of Fig. 
B). Therefore it appears to trace mostly the jet beam. Further out- 
wards, Peak B2 is centered in a hole behind the bowshock rim 
delineated in the near-IR image, while Peak B3 appears to lie on 
the bow rim (note however that the spatial extent of our SL map 
does not fully cover the tip of this region). 
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Fig. 3. (Ipper panel) H2 v=0-0 S(5) map (green contours) su- 
perimposed over the high velocity CO 7 = 2-1 map of Gueth 
|& Guilloteau| ( |1999] l; (Lower panels) Inner part of the H2 v=0- 
S(5) map (green contours) superimposed on the SiO J = 8-7 
and the CO / = 3-2 maps of |Lee et al. ( 2007| l. The grayscale 
backgrounds show the H2 v=l-0 S(l) image from McCaughrean 



et al. 



( |1994| >(top panel) or |Hirano et aL] ( |2006[ >(lower panels) 
The slight S-shape pattern of the S(5) emission map follows the 
brightness asymmetry of the cavity seen in the 2.12 pm image. 
In the blue-shifted jet lobe, point Bl is coincident with peak BII 
in CO 7=2-1 and BK4-6 in SiO. 



Symmetrically on the redshifted side, the mid-IR H2 emis- 
sion close to the source forms an extended curving "finger" trac- 
ing both the jet and the cavity wall north of it, as seen in the near- 
IR (see Fig. [3) and p eak Rl coincides with the near-IR knot F of 
|McCaughrean e t al. ( 1994). However, peak R2 is again centered 
on the cavity behind the bright bow rim traced in the near-IR. 
The spatial offset between B2 and R2 with respect to the near-IR 
counterparts suggests that mid-IR emission in these broad re- 
gions may be tracing the outer bow-shock wings where gas is 
expected to be in lower excitation conditions. 

The top row of Figure |4] presents the emission maps of the 
S(0) - S(l) rotational H2 lines observed with the LL - IRS mod- 
ules with a 10.5" sampling. Despite the lower angular resolution 
compared to the SL module, the HzS(Y) line clearly delineates 
the same outflow pattern as in Figure [3] As discussed in Sect. 
2, faint extended diffuse emission is also present, and becomes 
mostly evident at the north end of the map; diffuse emission is 
even more apparent in the H2 S(0) line map of Fig. |4] 

This diffuse brightness could be caused by the presence of an 
extended photo-dissociation region (PDR) created at the cloud 
surface by the illumination of FUV radiation from nearby stars. 
The existence of a diffuse PDR is strongly supported by the de- 
tection of diffuse PAH emission, which signals the presence of 
UV radiation in the field. The H2 diffuse emission has a bright- 
ness of the order of 10~ 5 erg/s/cm 2 for both the S(0) and S(l) 
lines: according to the Kaufman et al. ( 1999) model, such an in- 



tensity is compatible with a PDR with a density of 10 cm or 
larger and a FUV field of the order of 10 2 times the average in- 
terstellar field, measured in units of Go (i.e. in units of 1.6x10 3 
erg cirr 2 s~'). The most likely candidate for producing the dif- 



fuse PDR is the Per B0.5 star, located at a projected distance 
of about 15' north of the HH21 1 system. The FUV field of this 
star (assuming an effective temperature of 25 000 K, Hernandez 
et al. (2005)), diluted for the distance, can account for a Go ~ 
500, which is compatible with the estimation given above, tak- 
ing into considerations also projecting factors and partial field 
absorption along the path. We cannot exclude that part of the ex- 
tended S(0) and S(l) H2 emission also originates from other out- 
flows located north and south of HH21 1-mm, seen in both opti- 
cal and near-IR narrow band imaging of the region ( Walawender 
|et al.|2006{|Eisloffel et al.|2003b[ ). However, no close spatial cor- 
respondence is evident at our spatial resolution. 

3.2. Atomic emission 

The bottom row of Figure [4] presents the emission line maps of 
[SiII]35/im, [SI] 25jum, and [FeII]26jum obtained from the LL - 
IRS modules with a 10.5" sampling. The emission from atomic 
and ionic lines is substantially weaker in comparison to the H2 
lines. However, unlike H2 lines which are not detected close 
to the driving source, all three fine structure transition contours 
are continuous, evidencing emission very close to the protostar 
HH211-mm. Morphologically, the [SI] and [Fell] maps peak at 
the blue lobe bow-shoclQ Emission from [Sill] instead shows 
strong peaks on both the blue and red bow-shocks. 

The enhanced atomic - ionic emission on the blue-lobe bow- 
shock evidences a highly energetic interaction of the underlying 
jet with the ambient material. In the same region, a wealth of 
other molecular and ionic lines such as OH, H2O, HD, [Nell] 

2008] ) 



were mapped with Spitzer by Tappe et al 

rf 



suggesting a 

strong UV flux typical of a 40 km s~ l dissociative shock. The 
presence of such high velocity shocks in this region is also con- 
sistent with the detection of optical H a and [SII] emission spots 
(Walawender et al.[2006 2005). In the following, we will show 



that rotational H2 lines over the HH 211 outflow appear to trace 
slower shocks, either internal to the jet or in bowshock wings. 

4. Excitation conditions and mass fluxes 

4.1. H 2 emission 

4.1.1. Ortho to para ratio, extinction, and temperature 
components 

The H2 pure rotational transitions are easily thermalized in pro- 
tostellar outflow environments owing to their low critical den- 
sities (~10 2 -10 4 cirT 3 ) and can be used to probe temperatures 
in the range between ~300 and 1500 K. The most direct way 
to estimate the temperature of the emitting gas is by means of 
an excitation diagram (EXD); this involves plotting the quanti- 
ties ln(N u j/g u ) against E u j, where N u j and E u j are the column 
density and the energy of the upper level, and g u the statistical 
weight including the spin degeneracy (2S+1). Considering LTE 
conditions, the values of these quantities for the observed H2 
transitions should fall on a straight line, from the inverse slope 
of which the excitation temperature can be derived. Deviations 
from this linear correlation can be used to constrain other physi- 
cal properties, such as the H2 ortho-to-para ratio (hereafter, OTP) 
and the optical extinction Ay, as described below. 



1 In the [Fell] 26//m map, contours are interrupted just before the 
actual blue-shifted bowshock peak; this is an artifact resulting from a 
rogue pixel in the LL-IRS module that had to be masked out since it 
was introducing unreliable brightness values forming a bright vertical 
stripe across the map. 
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Fig. 4. Intensity maps of the HH211 flow in the PLSfl) and S(0) lines (top row) and the atomic [SiII]35ji/m, [SI] 25jum, and 
[FeII]26/im lines [bottom row) obtained with the Spitzer IRS - LL module (pixel scale ~10.5"). Crosses mark the position of 
the driving source HH21 1-mm. The dashed circle in the H2 maps masks the area around a bright star (IC348 IR) where line inten- 
sity measurements are unreliable. In the [Fell] 26pm map, contours are interrupted just before the eastern bowshock peak due to 
a rogue pixel. Contour levels in the H2 S(l) map start from 10~ 12 W cirT 2 sr~' with 10~ 13 increments, and from 7 10 ~ 13 W cm~ 2 
sr~' with 5 x 1CT 14 increments for the S(0) map. Contours of the ionic line maps start at 10~ 13 W cirT 2 sr~' with 10~ 13 W cirT 2 sr~' 
steps for [Sill] and [Fell], and at 1(T 13 W cirT 2 sr _1 with 2 x 1(T 13 W cirT 2 sr _I steps for [SI]. 



Deviations of the OTP from its LTE value are reflected as 
vertical displacements between the ortho and para transitions 
in the EXD, forming a "saw-tooth" pattern between the two 
H2 species. The observed OTP was estimated by examining the 
alignment of the S(5) data point with the neighboring S(4) and 
S(6) transitions, following the method outlined in Wilgenbus 



et al. (2000i; in all cases the spatial OTP variations are small, 



and values lie very close, within the statistical error limits, to the 
high-temperature LTE value of 3. 

The H2 S(3) transition at 9.7 pm is sensitive to the amount of 
dust along the line of sight, as it is located within a wide-band sil- 
icate absorption feature at the same wavelength; consequently a 
visual extinction value can be estimated by examining the align- 
ment of the S(3) point in comparison to the S(2) and S(4), having 
previously corrected for any deviations of the OTP ratio from the 
equilibrium value. Visual extinction is then estimated assuming 
an Agj/A v ratio equal to 0.087 ( |Rieke & Lebofsky|[T985) . Ay 



is found to exhibit no regular pattern along the outflow and to 
range between 7 and 9 mag, being consistent with previous es- 
timates of 7-15 mag by | McCaughrean et a l. (1994); Caratti o 
Garattietal.| ( |2006| r 



Considering an average value of Ay equal to 8 mag across 
the map, we have dereddened all the detected H2 lines falling in 
the SL range (ie. S(2) to S(7)) using the extinction law of Rieke 
|& Lebofsky ( |1985] l. Consequently, for all the points of the SL 
data cube where at least 3 H2 lines are detected above 3-sigma, 
we have calculated the excitation temperature from the slope of 
a least square fit on the dereddened data points in the EXD. An 



O'Connell et al. 



{2005} reported a higher extinction of 25 mag in the 



redshifted lobe but considered it rather tentative given their observing 
conditions 



example of these operations is shown in the left panel of Fig. [5] 
in the case of the Rl peak. 

The results over the whole HH21 1 region are summarized in 
Fig. [6] in the form of an excitation temperature map. Common 
morphological characteristics between the excitation tempera- 
ture map and the H2 S(2)-S(7) emissions are evident; peaks in 
excitation temperature tend to coincide with peaks of emission. 
At those peaks, temperature reaches values up to 1200 K, while 
in the more diffuse outflow regions the temperature is as low as 
700 K. In Table [T] we report the derived excitation temperatures 
with this method for the peaks of the S(5) emission along the 
blue and the red lobes of the outflow. 

In addition to the "warm" H2 component at ^ 1100 K indi- 
cated by the S(2)-S(7) lines, there is also evidence for a "cool" 
component from the S(0) and S(l) lines detected in the LL mod- 
ule. This is illustrated in the central panel of Figure |5J which 
present an EXD of the R1+R2 region in the LL 10.5" pixel, in- 
cluding the S(0) and S(l) datapoints. For consistency, the fluxes 
of S(2) to S(7) lines from the SL module were re-extracted at 
the coarser LL module sampling. Due to the extended diffuse 
emission, the S(0) line flux was obtained after subtracting the 
ambient contribution estimated in adjacent off-outflow regions. 

While the resampled S(3) - S(7) emission keeps the same 
"warm" temperature ~ HOOK, the data points for the S(0)-S(2) 
transitions show a steeper slope which is evidence for a "cool" 
component at T~ 300 K not probed by the higher J transitions. 
Another interesting result is that column densities obtained from 
the S(2)-S(7) lines in the LL pixel are smaller than in the SL 
pixel of 3.5", indicating that the warm component does not fill 
the LL pixel. 
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Fig. 5. Excitation diagrams of H2 in various pixel sizes. {Left) H2 columns in peak Rl from transitions in the SL module in 3.5" 
pixels (open circles). Ortho to para ratio and extinction are estimated from the alignment of the S(5) and S(3) data points with their 
neighboring S(2), S(4) and S(4), S(6) respectively. Dereddened columns (filled red circles) are linearly fitted (solid line) to derive 
the excitation temperature. Derived values are listed on the upper right corner of the diagram. {Center, right) Excitation diagrams of 
peaks R1+R2 and LL1 in 10.5" pixels, combining the H2 transitions observed in both the LL and SL modules. The SL line maps 
were re-extracted at the LL map pixel scale. A second cooler component is apparent in S(0) - S(2) with T\ ~ 300 K; the higher H2 
transitions trace a warmer component with temperature To ~ 1 100 K identical to that inferred at the SL pixel scale (see left panel). 



the data points in the EXD with the N/g axis, multiplied by the 
partition function Z{T) = 2 giexp{-Ej/kT). The later is esti- 
mated for the given temperatures, summing up for the first 35 
energy levels of H2. Derived column densities for the two tem- 
perature components are reported in Tables [T] and [2] and range 
from ~ 10 19 for the warm component (averaged over the SL 
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Fig. 6. Excitation temperature map of the "warm" H2 component 
deduced from the SL module map (3.5" sampling). Temperature 
is derived by a linear fit to the dereddened columns from the 
v=0-0 S(2) to S(7) lines; Peaks of S(5) emission from Fig.[2]are 
indicated. Pixels with less than 3 detected H2 transitions are in 
black. 



A similarly resampled EXD towards the LL1 peak also re- 
veals a cool component at 300 K, in addition to the warm com- 
ponent traced by the J > 2 lines, (see right-hand panel of Figure 
[5]). Such resampling of the SL data was not possible at the LL 
peak B2+B3, as part of the corresponding LL pixel falls beyond 
the SL map coverage. This double temperature structure differs 
from a similar Spitzer analysis of the L1448 outflow where all 
v=0-0 lines from S(0) to S(7) are interpreted with a dominant, 
single temperature component at 600-900 K. 

4.1.2. Column Densities and Mass Fluxes 

The total H2 column density of the "cool" and "warm" compo- 
nents can be measured from the intersection of the linear fit to 



pixel size) to 8 x 10 cm for the cool component (averaged 
over the LL pixel size). 

As already mentioned in the previous section, the inner peak 
of H2 emission at B1/LL1 coincides with peaks of the underlying 
CO and SiO jet; it is therefore instructive to compare H2 column 
densities with those deduced from CO, taking into account the 
differing beam sizes, to investigate possible physical connections 
between the different tracers. 

|Lee et aL] ([2007 ) estimate an H2 column towards the inner 
CO jet of 4 - 8xl0 z0 cirT 2 in their 1" beam, assuming an ex- 
citation temperature of 100 K and a standard CO/H2 abundance 
ratio of 8.5 x 10" 5 (ie. a fully molecular gas). If the H2 mid-IR 
emission arises from the same narrow jet, and is relatively uni- 
form along the jet axis, our derived column densities have to be 
corrected for the different beam sizes across the jet by a factor 
= pixel size / 1". This increases the H2 columns to ~ 3.5 x 10 19 
cirT 2 for the "warm" component, and to ~ 8 x 10 20 cm -2 for 
the "cool" component. The latter value is in excellent agree- 
ment with the CO-derived one. This coincidence suggests that 
the "cool" H2 gas detected by Spitzer towards LL1 could be trac- 
ing the same material as the CO jet, which would then be indeed 
mostly molecular. Note also that the temperature of 300 K of 
the cool H2 matches that inferred from SiO line ratios in the jet 
dHirano et al.(2006] >. 

In contrast, the warm H2 component column density cor- 
rected for beam dilution across the jet remains 10 times smaller 
than that of the CO jet. This may indicate that the warm H2 emis- 
sion close to the source arises only from a lower-density outer 
envelope around the CO jet. Alternatively, the warm H2 com- 
ponent could trace thin shock-heated zones within the CO jet, 
thus filling only a fraction of the 3.5" SL pixel length along the 
jet and suffering further beam dilution. We investigate the shock 
hypothesis in more detail in the next section {i 4.2 1. 

We may also determine a rough estimate of the one-sided 
mass flux of the cool and warm H2 components, assuming a uni- 
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Table 1. Physical properties of the warm H2 component from H2 v=0-0 S(2) - S(7) lines extracted at the SL (3.5") scale 



Position 


Offsets 


H 2 0-0 S(5) 


T 


N(H 2 ) 


M° 




(") 


(10- 11 Wcnr 2 sr" 1 ) 


(K) 


( 10'" cm- 2 ) 


(10- 7 M o yr-') 


Rl 


[-20.1, 13.3] 


3.3 


1100±70 


1.15 


1.1 


R2 


[-29.9, 17.0] 


2.3 


1050±75 


1.3 


1.2 


Bl 


[12.2, -4.7] 


1.3 


860±40 


1.1 


1.1 


B2 


[26.9, -10.3] 


2.0 


950±70 


1.5 


1.4 


B3 


[35.1,-10.6] 


1.9 


980±75 


1.4 


1.4 



M is proportional to (V^/lOO km s ')x(3.5"//,) and is not corrected for postshock compression. 



Table 2. Physical properties of the cool H2 component from H2 v=0-0 S(0) - S(2) lines extracted at the LL (10.5") scale 



Position 


Offsets 


H 2 0-0 S(l) 


T 


N(Ht) 


M" 




(") 


(10- 12 W cm" 2 sr" 1 ) 


(K) 


( 10 19 cm " 2 ) 


(10- 6 M o yr- 1 ) 


R1+R2 


[-19.3, 15.1] 


1.4 


330±50 


7.5 


2.0 


LL1 


[6.3, 2.0] 


1.3 


300±45 


7.8 


2.8 



" M is proportional to (Vj^/lOO km s ')x(10.5"//r) and is not corrected for postshock compression. 



form laminar flow across the corresponding pixel, using the re- 
lationship from Dionatos et al. ( 2009| l: 

M = 2fim H X{N(H 2 )A)x(V,/l t ) (1) 

where fi =1.4 is the mean weight per gram of hydrogen, mn is 
the proton mass, N(H2) the column density averaged over the 
area A of the pixel, /, is the projected emitting length along the 
flow direction (assumed in Tables 1,2 to be equal to the pixel 
length), and V, the tangential flow speed. 

The value of V, is quite uncertain. Taking into account radial 
velocities in the range 5-20 km s" 1 measured for the near-IR H2 



knots in the study of Salas et al. ( 2003[ ), and assuming an inclina- 
tion angle between 5° and 10° (see |TJ, tangential velocities may 
vary over 30 - 230 km s -1 , introducing an absolute uncertainty of 
a factor 3 either way in the inferred M values. In the following, 
we adopt V, = 100 km s _1 for ease of comparison with earlier 
work, bearing in mind that M values scale proportionally to the 
assumed velocity. In column 6 of Tables [T] and [2] we report the 
resulting M for the peaks of emission along the outflow, for both 
temperature components. 

For the cool H2 gas component, the laminar mass-flux is 
~ 20 - 28 x 10- 7 M Q yr -1 . For comparison, the one-sided jet 
mass- loss rate based on CO(3-2) emission obtained by |Lee et al.| 
(2007) assuming the same velocity of 100 km s -1 and a com- 
pression factor of 3 is 3.5 - 7 x 10 -7 M yr -1 . Adjusting their 
estimations for Helium, and uncorrecting for compression, the 
laminar mass flux of the CO jet rises to 15 - 30 x lO -7 M yr , 
this is in excellent agreement with the cool H2 laminar mass flux 
estimated here, again consistent with both tracing the same fully 
molecular gas. 

The warm H2 laminar mass-flux is 15-30 times lower at 
~ 1O- 7 M yr -1 x(V,/100km s" 1 ). Similar mass flux estimations 
were obtai ned from the pure ro tational H2 lines along the outflow 
of L 1448 dDionatos et al.|2009[ 0.5 - 1.4 X 1O- 7 M yr 1 ), where 



similar H2 columns were measurecj^j Note however that if the 
warm H2 does not trace a laminar flow as assumed, but arises in 
a single shock within the pixel, then the relevant emitting length 
l t along the flow will be generally narrower than the SL pixel 



size, and the true mass flux will be higher than listed in Table [T] 
Mass-fluxes in the warm component in the shock hypothesis are 
presented in the next section. 

4.2. Shock models for the warm H 2 component 



Note that a typo graphical error was introduced in Table 2 of 

are ~ 7 x 10 18 cm 2 at CS 



|Dionatos et al. (2009); the values of N(H 2 ) £ 
and ~ 1.5 — 3.8 x 10' 9 cnr 2 for the other pos 



positions. 



As pointed out in 54.1 peaks of warm (J>2) H2 emission coin- 
cide with peaks of temperature (Fig. [6]) and possibly correspond 
to regions of shocked gas. In order to constrain the shock condi- 
tions in these regions, we employ the existing shock model grid 
described in the work of Kris tensen et al.| ( [2007[ ) which is based 
on the MHD-VODE multi-fluid steady shock code of |Flower &| 
|Pineau des For ets (2003 ). The grid includes both continuous (C) 
and jump (J) type shocks and predicts the H2 lines brightness for 
various values of pre-shock density (n^'X shock velocity (V s ), 
initial OTP ratio (otp' m ) and transverse magnetic field density 
((p). The ranges of parameters investigated here are as follows: 

- n'™ = 10 4 to 10 7 cm - 3 by alternating factors of 5 and 2. 

- V s = 10-50 km s -1 with a step of 1 km s -1 

- otp mi = 0.01, 1, 2, 3 

- <p(f*G) - b x (cm -3 ) where we fix b- 1 for C and b-0 
for J-type shocks. 

The model is one-dimensional, considers 9 elements (H, He, 
C, O, N, S, Si, Mg and Fe) and 136 species connected with 1040 
reactions, and takes into account grain sputtering and erosion, as 
well as a variety of cooling and heating processes. Initial chem- 
ical abundances are computed assuming no UV field and an H2 
cosmic ray ionization rate of 5 x 10 -17 s _1 , yielding an ionization 
fraction ~ 7 x 10- 8 («™/10 4 cm s" 1 ) -0 ' 5 . The H 2 level populations 
and line brightnesses are integrated along the postshock cooling 
zone down to a temperature of 50 K, which is reached in a cool- 
ing time ~ 100 (10 s cm - 3 / «h) yrs in C-shocks (much shorter in 
J shocks). For more information on both the model and the grid 
data, the interested reader is referred to the articles of Flo wer &] 
|Pineau des Forets| ( |2003) and |Kristensen et al |d2007) ). 

In the following, we have employed a x l fitting method to 
optimally reproduce the observed H2 emission selecting the best 
matching C or J-type steady shock models from the grid. The x 1 
- fitting is performed with and without considering the S(2) tran- 
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Fig. 7. Steady ID C-shock model fits to the excitation diagram 
of the S(2)-S(7) lines at the peak of emission Rl in a 3.5" 
pixel. Observed points (filled circles with 10% error bars) are 
displayed along with best fitting x 2 models for pre-shock den- 
sities between 10 4 and 10 7 cirT 3 . The top panel lists the corre- 
sponding shock velocities and beam filling factors, along with 
the reduced x 2 fit values to the S(2) - S(7) (J>1) and to the S(3) 
- S(7) (J> 2) H2 transitions. A strong degeneracy among models 
is found (see text). 



sition, since the latter may be dominated by the cool component 
(see Fig|5]). 

4.2.1. C-Shocks 

As a general trend we found that steady C-shock models are able 
to reproduce the observed S(2)-S(7) brightnesses only if we as- 
sume a small beam filling factor. For this reason we have per- 
formed our^ 2 fits on relative line brightnesses, scaled to the H2 
S(5) line; for the best fitting models we have then calculated 
a beam filling factor as the average ratio of the observed over 
model line brightness. 

In Figure[7]we plot the EXDs of the best fitting;^ 2 models for 
each preshock density, scaled by the corresponding beam filling 
factor, along with the observed data (filled circles) for the peak 
of emission Rl; note that beam filling factors decrease for in- 
creasing values of density. There is a strong degeneracy among 
the models, with a very small spread in column densities for n™ 
in the range between 5 10 s and 10 7 crrT 3 and velocities in the 
range 10-13 km s _1 , however models for lower values of n# < 
10 5 and velocities up to 30 km s _1 cannot be excluded; each one 
of these models can account fairly well for the observed emis- 
sions. 



Figure [8] presents EXD of the best matching C-type shocks 
(red squares, always scaled by the corresponding beam filling 
factor) along with the observed data points (filled circles), for 
the peaks of emission as observed along the outflow. Input pa- 
rameters for models are reported in Table|3] All C-shock models 
predict moderate shock velocities (10-15 km s _1 ) and an initial 
ortho to para ratio at the high temperature LTE value of 3. 

Steady C-shock models reproduce the data only if we con- 
sider low beam filling factors bff ~ 0.01-0.04, or in other words 
a shock area A s h oc k -A p i xe iX bff, corresponding to a typical size 
of ^ A s h ck ~03" - 0.7". These dimensions are comparable to 
the reported jet width <\" in the interferometric CO and SiO 



maps of |Lee et aL ( |200 7). High post-shock densities in the same 



range as the best fitting C-shock models (5 10 5 - 10 7 cm 3 ) have 



been retrieved in the jet by SiO emission studies fNisini et al. 
2002b{ pirano et al.|2006||Lee et al.|2007[ ), supporting our mod- 
eling results. The shock models predict thin H2 post-shock cool- 
ing zones of 0.015"- 0.25", much smaller than the SL pixel size 
of 3.5", explaining the small warm column densities derived in 
the previous section. The cooling zone is also smaller than the 
shock surface size inferred from the filling factor, as required for 
the 1-D planar shock approximation to be locally valid. 

If the warm H2 emission indeed traces low-velocity internal 
shocks within a fast jet beam, we may infer the jet mass-flux 
from the pre-shock density and the shock area as: 



M j<shock = fiiriH x bffApbei x ri^" x V jet . 



(2) 



with a beam 

-1 



Considering a pre-shock density of 5 x 10 cm 
filling factor bff of 0.023, and a jet speed of 100 km s~\ we 
obtain for the warm H2 Mj tS hock = 7 x 10~ 7 M Q yr _1 . This is 7 
times higher than the value obtained in Section 4.1.2 assuming 
a laminar flow filling the SL pixel. The mass-flux entering the 
shock is only a factor 3-4 below the laminar mass-flux for the 
cool H2 (and CO). Such a discrepancy is expected in a time- 
variable flow undergoing internal shocks (Harti gan et al.|1994 i. 
Therefore, the "warm" H2 component seems compatible with 
localized heating in slow internal C-shocks within the molecular 
jet beam, and may not trace a separate lower density flow. Note 
that pre-shock densities higher than 5 x 10 6 cirT 3 seem to be 
excluded, as the pre-shock mass-flux rate from Equ. [2] would 
become comparable to the accret ion rate onto th e HH211-mm 



Lee et al. 2007 



driving source (8 x 10 6 M Q yr~ 

However, SiO is observed only within ~ 15" from the source, 
therefore pre-shock densities of 5 x 10 5 cm 3 might apply only to 
the Bl region. In addition, in our maps the H2 emission remains 
unresolved only close to the driving source. Further out at the 
bow-shocks, the emission appears slightly resolved at the 3.5" 
scale. A small bff < 0.04 would therefore imply that the bow- 
shock surface is highly patchy. This is a reasonable possibility, 
since HST observations in several outflows have shown that what 
appears as an individual bow-shock at low resolution is in fact 
resolved in many sub-structures of sub-arcsec scale (e.g. HH2 in 
Orion, |BaIiy et al.| ( |2002] >). In H H211, SMA obse rvations, per- 
formed with a resolution of ~ 1" ( |Lee et al.|2007| l, show indeed 
that the CO emission is highly clumped, both in the high velocity 
jet and along the low velocity cavity ahead of the bow-shocks. 

An alternative explanation for the small beam filling factor in 
the outer HH21 1 bow-shocks might be that the pre-shock density 
drops to < 10 4 cm -3 , so that the cooling length becomes compa- 
rable to the 3.5" pixel size; the hypothesis underlining our model 
comparison, that the whole cooling zone fits inside the pixel, 
would then become invalid. Indeed, [O'Connell et al.| ( |2~005| l de- 
rived a low pre-shock density in the Rl and R2 regions by fitting 
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Fig. 8. Excitation diagrams for all emission peaks mapped with the SL module; observed points (filled circles with error bars), are 
displayed along with best^ 2 fitting C (red squares) and J (blue triangles) shock models. C- shock models can account fairly well 
for the observed emissions considering a small beam filling factor, whereas J-shock models cannot reproduce well the slopes of the 
observed data points, except marginally in Rl and R2. 



a 3D C-type bow-shocks to the near-IR ro-vibrational H2 lines. 
However, no predictions for the mid-IR H2 emissions were pro- 
vided. Extension to 3D bow-shock geometry of the model ap- 
plied here (Gustafsson et al. 2010) > will be needed to perform 
more reliable modeling of the outer bowshocks of HH21 1 than 
possible in the present analysis. 

A third possibility might be that the preshock gas ahead of 
the leading bowshocks is more ionized than we have assumed, 
due to the presence of a UV flux from nearby dissociative shocks 
such as that observed by Tappe et aT] ( |2008| ) in the blue-shifted 
apex. This would increase the ion-coupling in the C-shock, re- 
ducing its width and column density, and yielding a larger filling 
factor. Indeed, [O'Connell et al.| ( [20 05 ) find that a high pre-shock 
ionization fraction of 10" 3 and a high atomic hydrogen fraction 
of 20% provide a good fit to the near-IR lines in the redshifted 
bow. However, given the already high degeneracy of our model 
fits, and the lack of information on shock geometry and thick- 
ness in low angular resolution Spitzer maps, introducing the UV 
field as an additional free parameter is not feasible in practice. 
Additional constraints from other line tracers or from higher res- 
olution mid-IR maps (eg. with JWST) will be required to further 
explore this alternative. 

4.2.2. J-Shocks and non-steady C-shocks 

In the case of J-shocks, absolute brightness can be fairly well 
matched by models and beam filling factors < 1 were not nec- 
essary to be considered. The predicted pre-shock densities for 



J-shock models ~ 1 - 5 X 10 4 crrT 3 are lower than for the best 
fitting C-shock ones. As a general trend though, the best fitting 
J-shock models, presented with blue triangles in Fig.[8]and listed 
in Table|4] cannot reproduce at the same time the slope of the ob- 
served data points, except for the red-shifted bowshock (Rl and 
R2). 

An analogous behavior has been noted in the work of 
Giannini et al.|(|2006|l, where mod els produced by the same code 
of Flower & Pineau des Forets ( 2003 \ are employed to study 
the conditions of the HH54 outflow. Considering only the ro- 



tational H2 emission in reported excitation diagrams, Giannini 



et al. ( 2006 1 find that steady C-type shocks overpredict the ob- 



served v=0 column densities, whereas J-type shocks can fall 
much closer but cannot fit at the same time the observed slope. 
In the case of HH54, it is concluded that a non-steady, truncated 
C-shock with an embedded J-type front (C+J shock) optimally 
reproduces the observed H2 emission. The same conclusion was 
reached also in the case of LI 157 ( |Gusdorf et al.|2008a| ). Indeed, 
the overall effect as reflected in excitation diagrams is to lower 
the column densities of low excitation energy levels (truncated 
C-precursor) and increase the column densities of the rovibra- 
tional energy levels of H2, excited in the J-type front ( |Flower| 
et al.|2003] l. In our case, such a model might potentially fit our 
mid-IR data in the outer bowshocks without the need to con- 
sider a beam filling factor < 1 . A more detailed comparison with 
unsteady C+J-shock models would be needed to test this pos- 
sibility, however such models were not included in our adopted 
grid as the age of the shock would have, again, added an extra 
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Table 3. Parameters of shock models best fitting the S(3) to S(7) lines 



Peak 


n™(cnT 3 ) 


C-shock (b= 
V s (kms-') 


= D 
opr,„, 


b ff 


J-shock (b = 0, bff= 
n™(cnT 3 ) V^kms- 1 ) 


1) 

opr,„, 


Rl 


5 10 5 - 10 7 


10- 13 


3 


0.01 - 0.02 


5 10 4 


14 


3 


R2 


5 10 5 - 10 7 


10-11 


3 


0.02 - 0.04 


5 10 4 


14 


3 


Bl 


5 10 5 - 10 7 


10 


3 


0.02 - 0.03 


10 4 


10 


3 


B2 


5 10 5 - 10 7 


10-11 


3 


0.01 - 0.03 


10 4 


10 


3 


B3 


5 10 s - 10 7 


10- 11 


3 


0.02 - 0.03 


10 4 


11 


3 



free parameter to already degenerate fits. The shock ages would 
have to be much shorter than the outflow dynamical age of about 
1000 yrs in order to reduce by a factor 50 the warm H2 column 
(typically a shock age of 1/50 of the H2 cooling time, ie. ~ 2 
(10 5 cm ~ 3 / n H ) yrs). 



4.3. Atomic emission 

As pointed out in section ^3] all atomic/ionic emission falls 
within the LL IRS module and consequently these maps are of 
inferior spatial resolution compared to the SL ones. As atomic 
emission is found to be enhanced at the bow shocks but remains 
significant also close to the driving source, we have focused our 
analysis all along the jet axis including the intermediate points 
LL1 and LL2 (see Fig. [4} between the bow shocks, where the 
[Fell] 26/im line is detected. For the analysis of the ionic emis- 



sion we have followed the techniques employed by Dionatos 
|et al.| ( |2009| l, and the reader is referred to this article for a de- 
tailed description. 



0.8 
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4.3.1. Electron density 

As a first step, the observed emissions of [Sill] and [Fell] are em- 
ployed in order to constrain the electron density. In Figure [9 the 
ratio of [SiII]34.8yum over [FeII]26.0jt/m as a function of the elec- 
tron density is presented, for temperatures of 500K, 1000K and 
4000K (solid and dashed lines). A statistical equilibrium model 
was employed that uses the first 16 levels for [Fell] (Nisin i et al.| 
2002a), and a two level system was considered for [Sill], with 
radia tive and collisional rates with electrons taken from lDuftonl 
& Kingston ( 1991[); a solar abundance ratio of Si + /Fe + was con- 



sidered, taking into account the similar ionization potential of 
the two species as well as the similar efficiency of sputtering of 
Si and Fe atoms from olivine grains in C-shocks, according to 
Figure 2 of Gusdo rTet al.| ( |2008b) . In the diagram we plot the 
observed ratio limits at the bow shocks for clarity, but the re- 
sults for the other points along the jet are similar and are listed 
in Table [4] Lower limits for the electron density vary between 40 
and 150 cirT 3 for temperatures of 4000K and 500K respectively, 
while upper limits converge for the considered temperatures to 
~ 400 cm 3 . Temperature limits lower than a few thousand K are 
reported in the work of Tappe et al. ( 2008) for the blue-shifted 
bow shock considering a non-LTE analysis of the three [Fell] 
lines (18 /urn, 26/im and 35yum) that are detected with the high 
resolution IRS modules. Based on this, we restrict the electron 
density lower limit to 100 cirT 3 by assuming a maximum tem- 
perature of 3000 K. 



Fig. 9. Diagnostic diagram of the [SiII]35jum/[FeII26]^m ratio 
versus electron density, for temperatures of 500K, 1000K and 
4000K (solid and dashed lines) and assuming a solar Si + /Fe + 
abundance ratio. The hatched area represent the observed ratio 
for the north and south bow-shock regions, from which electron 
densities between ~ 100 and 400 cirT 3 are found . 



4.3.2. Gas phase depletion 

We investigate the depletion onto grains following the analysis 
of Dionatos et al. ([2009). We calculate the gas phase abundance 
of refractory elements such as iron and silicon by comparing 
their emissions against the emission of a non-refractory species, 
in the present case sulphur. Theoretical emissivities for sulphur 
were calculated considering a five level statistical equilibrium 
code for an upper temperature limit of 3000 K assumed from 
the work of (Tappe et al. 2008) and a lower limit of 700 K as 
imposed from the H2 analysis. For the excitation of sulphur, we 
have also calculated emissivities due to collisions with atomic 
hydrogen, taking into account that in environments with low ion- 
ization fraction and high total density such excitation may prove 
to be significant, despite the fact that collisional de-excitation 
rates for atomic hydrogen are lower than the electronic ones 
( jHollenbach & McKee|[T989l y H ly e ~ 0.02) for the [SI]25jum 
line. For these calculations we have considered a medium with 
n(H) = 10 4 cm ~ 3 . Excitation of [Fell] and [Sill] by H-collisions 
was not taken into account because neutral hydrogen collisional 
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Table 4. Physical properties of the atomic/ionic line component 



11 



Position 


Offsets 
(") 


[SI] 

(10"' 3 Watt cm- 2 sr" 1 ) 


[FeII]/[SI] 


[SiII]/[FeII] 


M([SI]) fl - fc 
(lO^Moyr 1 ) 


b,c 
H 

(10 s cm" 3 ) 


R1+R2 


[-19.3, 15.1] 


5.34 ± 0.13 


0.92 ± 0.04 


2.18 ± 0.51 


0.4 - 3.4 


0.21 - 1.61 


LL1 


[-6.5, 8.6] 


6.24 ± 0.16 


0.60 ± 0.02 


1.77 ±0.82 


0.4-4.1 


0.24-1.88 


LL2 


[6.3, 2.0] 


7.49 + 0.19 


0.35 ± 0.01 


1.95 ± 0.93 


0.7-4.9 


0.29 - 2.91 


B2+B3 


[31.9,-11.0] 


16.1 ± 0.41 


0.37 ± 0.01 


2.14 ± 0.42 


1.3 - 10.5 


0.62 - 4.9 



"Mis proportional to (Vj^/lOO km s"')x(10.5"//r) and is not corrected for postshock compression. 
* values calculated for collisions with electrons at T = 3000 K (lower value) and T = 700 K (upper value); including collisions with hydrogen for 
n(H)=10 4 cirT 3 would decrease the minimum and maximum values by a factor of ~ 2 and 4 respectively. 
r average proton density assuming an emitting volume of 10.5"x l"x 1" 



Table 5. Gas phase abundance of Fe + and Si + from line ratios to 
sulfur 



1.5 



S l 
a. 

a; 
m 




0.5 



3. 

CD 

a; 



E, o 



-0.5 



-| — i — i — i — i — | — i — i — i — i — | — i — i — i — 

n B =100 cm- 3 

n =400 crrr 3 



_n e =400 cm" 3 , n(H)=10* cm" 



R1+R2 



B2+B3 



1000 2000 3000 4000 5000 
Temperature (K) 



Fig. 10. Diagnostic diagram of [FeII]26/mi/[SI]25/im ratio ver- 
sus temperature assuming solar relative abundances. Electron 
densities of 100 and 400 cirT 3 are considered (full and dot- 
ted lines, respectively). The effect of collisions with atomic hy- 
drogen on the [SI] excitation is also shown (long-dashed line, 
n(H) = 10 4 cirT 3 ). The observed ratios are always below the 
predicted values for solar abundances, which is interpreted as 
iron being depleted onto dust-grains in the considered regions. 



rates are much less important for ions (i.e. ~ 10 3 -10 4 weaker 
than electronic collisional rates), and are not available for all the 
levels considered in our model. 



Results of these calculations are illustrated in Fig. 10 where 



we plot the ratio of [Fell] 26yum over [SI] 25/itn against gas tem- 
perature for electron densities of 100 and 400 cm" 3 , assuming 
that iron is fully ionised (which is the case in most dissociative 
shocks) and has a solar abundance. The observed ratios for the 
points R1+R2 and B1+B2 are way below the expected values, 
giving rise to the hypothesis that iron is heavily depleted onto 
dust grains in the considered regions, more so if the emitting gas 
is cooler. The case including atomic hydrogen collisions with 
n(H) = 10 4 cm 3 is separately plotted in Fig. 10 The iron de- 
pletion is less severe in this case but is still substantial, about a 
factor 10. 



Position: 


R1+R2 


LL1 


LL2 


B2+B3 


[Fe s + a5 ]/[Fe G ]« 
[Pel, 1/lFea]* 


4-9 10" 2 
10-20 10" 2 


3-6 10" 2 
7-10 10" 2 


2-4 10" 2 
5-7 10" 2 


2-4 10" 2 
5-7 10" 2 


[si;j/[si r 
[si^j/tsy* 


4-9 10" 2 
10-20 10" 2 


2-4 10" 2 
5-8 10" 2 


1.5-3 10" 2 
3-6 10" 2 


2-4 10" 2 
4-7 10" 2 



" assuming only collisions with electrons (n e = 100 - 400 cm" 3 ) and 

r f =700 K (lower value) and 3000 K (upper value) 
b including excitation of sulfur by atomic hydrogen with n(H) = 10 4 



In Table [5] we give the inferred gas phase depletion for tem- 
peratures of 700 and 3000 K, both without and with atomic hy- 
drogen excitation of [SI]. The gas phase abundance of Fe + is 
about 2-20% of the solar value for T= 3000 and 700 K respec- 
tively, for all the points under examination. Similar low values 
were found also in the case of the L1448 jet (5-20%). Identical 
gas-phase depletions are obtained for Si + , but since we assumed 
a solar Fe + /Si + ratio in the first place to derive our n e values, this 
is just a self-consistency check. Interestingly, the SiO abundance 
in the HH21 1 molecular jet, estimated as ~ 10" 6 by (Nisini et al. 



2002b] |Chandler & R icher 2001) corresponds to 3% of the solar 
abundance of silicon, therefore the fraction of silicon released 
in the gas phase seems comparable in the atomic and molecular 
components of HH2 1 1 . 

4.3.3. Density and mass flux of the atomic component 

The total number of H nuclei in the emitting volume V can be 
derived from the line luminosity according to the relation: 



(X 1\ 
hvAifi — J , 



(3) 



where A, and /; are the spontaneous radiative decay rate and frac- 
tional population of the upper level, and [X/H] is the gas phase 
abundance of the atom/ion under consideration with respect to H 
nuclei. Line luminosities were integrated over the LL pixel size 
area, and the fraction of atoms at the upper level was calcu- 
lated for n e between 100 and 400 cm 3 and T= 700-3000K as 
imposed from the previous analysis. 

This relation was applied to [SI], the only non-refractory 
species detected, assuming that all sulfur is in neutral form and 
employing its solar abundance (Asplu nd et al.|20 05). Such cal- 
culations for refractory species would be redundant, as their gas 
phase abundances from Table [5] were determined from their in- 
tensity ratio to [SI]. 
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Values of rin at individual positions are listed in column 7 of 
Table [4] assuming an emitting volume in the LL pixel of 10.5"x 
l"x 1" (uniform narrow jet). The inferred n# depends strongly 
on the adopted excitation conditions: it is about 0.2 - 0.6 x 10 5 
cirT 3 for T=3000 K, n e = 100 crrT 3 , and 8 times higher for 
T=700 K, n e = 400 cm" 3 . The ionisation fraction is then 1.6 - 
5 x 10" 3 at 3000 K, or twice smaller at 700 K. Such an ionization 
level seems to favor the upper temperature, and therefore the 
lower density range. If the fraction of H atoms is important with 
n(H) a 10 4 cm" 3 , the values would be further reduced by 
a factor 2, down to a few 10 4 cm" 3 . In the case again, that the 
emission comes from a smaller volume, eg. a small 1" knot as 
observed in the near-IR [Fell] 1 .64/mi line, the density would rise 
up to 10 5 cm" 3 . 

The mass flux in the atomic/ionic component assuming a 
laminar flow can also be derived applying the relationship given 
inp 



i Nisini et al. (2005) as applied to the mid-IR lines in Dionatos 
|et al.| < |2(J09) : 



M = nm H x (n H V) x (V,/l t ), 



(4) 



where fi— 1.4 is the mean weight per H nucleus, m# is the proton 
mass, Yin V is the total number of protons in the emitting region 
as given by Eq. [3] I, is the projected emitting length along the 
flow and V, the tangential flow velocity. 

Derived mass flux values at the various emission peaks are 
listed in column 6 of Table [4] for a tangential speed of 100 km 
s" 1 and I, = 10.5" (LL pixel size). Like rinV, they again depend 
strongly on the adopted excitation conditions. From the ioniza- 
tion fraction considerations made above, we favor the high tem- 
perature, smaller M values of ~ 0.4 - 1 x 10" 7 M yr" 1 . These 
are 20-50 times smaller than the cool H2 jet mass-flux, assuming 
the same flow speed of 100 km s" 1 . 

However, we stress again that M estimates assuming a lam- 
inar jet flow may be in error if the [SI] emission arises from 
shocks, which is likely given the ionization fraction inferred 
above. This is difficult to quantify without appropriate shock 
models. 

The atomic/ionic emission appears to require higher excita- 
tion shocks than those producing the mid-IR H2 lines. The low 
velocity 10-15 km s" 1 shock models that best fit the warm H2 
emission are unable to reproduce also the observed [Fell] and 
[Sill] line intensities unless the shocks are of J-type and these 
atoms are essentially undepleted. As discussed in Section 4.3.2 
there is overwhelming evidence of depletion of these refractory 
species, so [Fell] and [Sill] probably arise from faster shocks 
than those dominating the mid-IR H2 lines. It is indeed likely 
that the mid-IR atomic lines originates from the same dissocia- 
tive shocks at v s > 30 km s" 1 that give rise to the optical and NIR 
emission of H g , [SII] and [Fell] (jO^Connell et al.|2005||Caratti| 
|o Garatti et al.|2006||Walawender et al.|2006||2005) . ~ 

Unfortunately, our shock models do not yet include ioniza- 
tion and dissociation by the shock UV flux, so we cannot explore 
the relevant excitation range. However, we note that the J-shock 
models of |Hollenbach & McKeel ( |1989| l, which include UV flux, 
suggest that the [SI]25/im brightness is relatively independent of 
shock speed over the range 30-100 km s" 1 and is roughly pro- 
portional to n™. Our measured flux inside the LL 10.5" pixel 
would suggest n™ x bff 100 - 1000 cm" 3 . From Eq. |5]we 
would infer a preshock mass-flux of ~ 0.5 - 5 x 10" 7 M Q yr" 1 , 
still smaller by a factor 4-40 than the cool H2 jet mass-flux. 

These results would suggest that the atomic component in 
the HH21 1 jet does not have enough momentum flux to entrain 
the CO/SiO/H2 jet, their velocities are comparable and the 



shock compression in the cool H2/CO component does not ex- 
ceed a factor 3. The molecular jet would then have to trace mate- 
rial ejected from the accretion disk, while the atomic component 
would trace a separate ejection, e.g. from hotter more internal 
regions of the accretion disk. On the other hand, a higher com- 
pression factor in the cool H2/CO emitting zone, or a faster ionic 
jet would suffice to remove the discrepancy, therefore a definite 
conclusion on whether the molecular jet is ejected or entrained 
cannot be reached from the present data alone. 

5. Conclusions 

We have carried out Spitzer spectral mapping observations 
towards the jet driven by the Class source HH211-mm. 
Molecular lines (pure rotational H2) as well as fundamental 
atomic and ionic lines ([SI], [Sill], [Fell]) were detected, and 
their maps follow the characteristic bipolar outflow pattern as 
traced by near-IR H2 and CO lines. H2 emission becomes im- 
portant only 5" away from the driving source while atomic and 
ionic lines are detected very close to the driving source. In the 
inner part of the blue-shifted lobe, the H2 emission is spatially 
coincident with the high velocity jet observed in CO and SiO. 

Analysis of the observed H2 lines reveal two temperature 
components: "cool" gas at T~ 300K dominating the mass, and 
"warm" gas at T~ 700-1200K of 10 times smaller average col- 
umn density. Once corrected for beam dilution, the column den- 
sity of "cool" H2 towards the CO jet is compatible with hydro- 
gen being mostly molecular. The warmer component traced by 
the S(2) to S(7) lines is well fitted by C-shock models of high 
density 5 x 10 5 and a small shock cross section of 0.5", com- 



patible with the density and width of the CO jet quoted by Lee 
et al. (2007). The warm H2 emission could then trace thin layers 
of warm post-shock gas within the CO jet. 

Similarly, high-density C-shock models can also account for 
the brightness of the H2 mid-IR emission further downwind be- 
yond the CO and SiO jet, but the small shock surface filling fac- 
tor of 0.01-0.04 is not easy to reconcile with the spatial exten- 
sion visible in Spitzer maps, unless the bow-shock wings are 
very clumpy. Lower density shocks of n# < 10 4 cm" 3 where the 
cooling zone spreads over several pixels and/or higher preshock 
ionization may need to be considered in these regions, as previ- 
ously invoked to model near-IR H2 excitation in the red-shifted 
bow-shock of HH21 1 ( P'Connell et al.|2005) . 

The detected fine structure lines mapped very close to the 
driving source signify the presence of an embedded atomic jet. 
Line ratio diagnostics indicate a gas-phase depletion of iron and 
silicon of at least a factor 10, and lower excitation conditions 
than in optically visible jets. An atomic jet of similar proper- 
ties has been also detected by Spitzer in the outflow of L1448-C 
( Dionatos et al.|[2009| >. As suggested in that case, the detected 
atomic gas in the HH21 1-mm outflow may represent the equiv- 
alent for the Forbidden Emission Line (FEL) region observed in 
more evolved ClassI/II sources. The gas depletion of iron and sil- 
icon indicates that dust grains have survived in the atomic flow, 
ruling out an origin from within the dust sublimation zone close 
to the protostar. The excitation conditions require faster shocks 
than those producing the H2 mid-IR lines. 

Estimations of the molecular and atomic mass-flux rates 
have been performed using both a laminar flow assumption and 
shock models. The cool H2 mass-flux is comparable to that in- 
ferred from CO observations, with a value uncorrected for com- 
pression of ~ 2x 10" 6 Mq yr" 1 (V/lOOkm s" 1 ); a similar value is 
found for the warm H2 if it arises from dense C-shocks as sug- 
gested by best fit models. The atomic component mass-flux is 
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uncertain by a factor 8 due to its uncertain temperature, but con- 
siderations of ionization fraction as well as published dissocia- 
tive shock models suggest a mass-flux smaller than in the molec- 
ular jet. However, the momentum fluxes could become compa- 
rable if the atomic jet is faster than the molecular jet, or if the 
cool H2 and CO suffered shock compression. Given the uncer- 
tainties, it is still unclear whether the molecular CO jet is tracing 
ambient gas entrained by the atomic jet, or if it has to be ejected, 
for example in a molecular MHD disk wind as recently modeled 
by Panoglou et al. (2010, A&A, submitted). 
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